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Abstract 

Using proton-proton collision data collected by the LHCb experiment at y/s = 7 TeV, 
corresponding to an integrated luminosity of 1.0 fb _1 , the ratio of branching fractions 
of the B° — > D*~7r + 7r~7r + decay relative to the B° — > D*~tt + decay is measured to 
be 

B(B° -> L>*-vr+7r-vr+) , , 

B(B°^D*-n+) = 2M ± ° M (8tat) ± °' 13 (sySt ° • 

The Cabibbo-suppressed decay B° — > D*~ K + ir~7r + is observed for the first time 
and the measured ratio of branching fractions is 

B(B° — > D*~ K + ir~ir + ) 

B \ B o r^Sr+^r+y = (6 - 4? ± °' 37 ^ ± °' 35 (sySt0) X ^ ' 

A search for orbital excitations of charm mesons contributing to the 
B° — > D*~ir + ir~ir + final state is also performed, and the first observation of the 
B° — > Di(2420)°7r + 7r _ decay is reported with the ratio of branching fractions 

B(B° -> (^1(2420)° ->• D*-vr+)vr-7r+) . , 9 

-i b{B0^D*-^-k+) = (2 ' 04 ± °- 42 (stat ° ± °- 22 (sySt - )} X 10 ' 

where the numerator represents a product of the branching fractions 
B{B° ^i(2420)°^-^+) and i3( J Di(2420)° -> D*-tt + ). 
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1 Introduction 



Open charm decays of b hadrons offer a means by which both the electroweak and QCD 
sectors of the Standard Model (SM) may be tested. Beyond measurements of CP violation 
and the phases derived from the CKM matrix, rare B — > DX decays may be used to 
search for new physics in decays mediated via annihilation or exchange processes. High 
multiplicity B — > DX decays are receiving increasing attention at LHCb j^j3], in part 
owing to the large samples that can be obtained as a result of the copious bb production 
at the LHC g|. 

Improving knowledge of the B° — > D^*^n + n^n + decay is of interest because of its 
potential use as a normalisation mode for the semileptonic decay B° — > D^*^t + v t with 
r + — > 7t + 7t~7t + v t |5|. Charge conjugation is implied throughout. The latter B decay has 
recently shown an excess over the SM branching fraction expectation [HJ, which could 
indicate the presence of physics beyond the SM. 

In this work three measurements of ratios of branching fractions are described. The 
first measurement is of the Cabibbo-favoured (CF) ratio 

B(B° D*-tt + it-it + ) 

r 3 h = 



B(B° D*-n+) 



from which a value for B(B° — > D* tt + tt 7r + ) is obtained. The current world average 
value is B(B° D*-n+TT-n+) = (7.0 ± 0.8) xl0~ 3 Q. The subscript 3h denotes the 
three hadrons (or "bachelors" ) in the signal B° decay, which are produced along with the 
charmed meson. The second measurement is of the ratio of branching fractions of the 
Cabibbo-suppressed (CS) decay B° — > D*~ K + n^n + relative to its CF counterpart 



Kirn 



r cs 



B(B° ->• D*-K + n-n + ) 
B(B° -> D*-ir+ir-n+) 



where the subscript denotes Cabibbo suppression of the signal decay. The third measure- 
ment is of the ratio of branching fractions of the CS decay B° — > D*~K + relative to its 
CF counterpart 

K _ B(B° -j- D*-K + ) 
r ° s ~ B{B° D*-ti+) ' 

The superscripts in r^ s and are used to differentiate the single- and triple-bachelor 
measurements, respectively. The rate of the CS processes are expected to be smaller than 
their CF counterparts by a factor of roughly tan 2 (6 l c) ~ 1/20, where 8q is the Cabibbo 
angle [fj. 



2 Detector, signal selection and simulation 

The LHCb detector [8] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r\ < 5, designed for the study of particles containing b or c quarks. The detector 



1 



includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream. The combined tracking system has a 
momentum resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 100 GeV/c, and an 
impact parameter (IP) resolution of 20 [im for tracks with high transverse momentum (pr)- 
Charged hadrons are identified using two ring-imaging Cherenkov (RICH) detectors, which 
provide a kaon identification efficiency of ~95% for a pion fake rate of a few percent, 
integrated over the momentum range from 3-100 GeV/c M. Photon, electron and hadron 
candidates are identified by a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons 
are identified by a system composed of alternating layers of iron and multiwire proportional 
chambers. 

LHCb uses a two-level trigger system. The first level of the trigger consists of a 
hardware stage that searches for either a large transverse energy cluster (E^ > 3.6 GeV) 
in the calorimeters, or a single high p^ muon or di-muon pair in the muon stations. This is 
followed by a software stage that applies a full event reconstruction. The software trigger 
requires a two-, three- or four-track secondary vertex with a high sum of the px of the 
tracks and a significant displacement from the primary pp interaction vertices. At least one 
track should have pr > 1-7 GeV/ c and IP \ 2 with respect to the primary interaction greater 
than 16. The IP \ 2 is defined as the difference between the x 2 °f the primary vertex (PV) 
reconstructed with and without the considered track. A multivariate algorithm is used for 



the identification of secondary vertices consistent with the decay of a b hadron 10 . The 
results presented in this paper use the full 2011 dataset corresponding to an integrated 
luminosity of 1.0 fb -1 of pp collisions at = 7 TeV. 

Candidate B° decays are formed by combining a D*~ — > (D° — > K + 7t~)tt~ candidate 
with either a single bachelor or triple-bachelor combination. All final state tracks are 
required to have in excess of 500 MeV/c 2 , except for the slow pion produced in the 
D*~ decay, which must have px in excess of 100 MeV/c. All tracks must be well separated 
from any reconstructed PV in the event. They must be identified as either a pion or a 
kaon using information from the RICH detectors. Particle identification likelihoods for 
several hypotheses (e.g. ir + ,K + ,p) are formed, and the difference in the logarithms of 
these likelihoods, ALL, is used to differentiate between particle types. Candidate D° 
mesons are required to have good vertex fit quality, be well displaced from the nearest PV 
and have an invariant mass m(K + 7T~) within 50 MeV/c 2 of the D° mass. Candidate D*~ 
decays are selected by requiring 140 < m(K + 7i^7i + ) — m(K + 7r~) < 150 MeV/c 2 . 

Candidate triple-bachelors are formed from a tt + 7t~tt + or K + tt~tt + combination, where 
all invariant mass values up to 3 GeV/c 2 are accepted. The vertex of the combination 
must be well separated from the nearest PV. Backgrounds from B° — > D*~7t + (tt~7t + ) 
decays for the CS modes B° — > D*~ K + (-k~ -k + ) are reduced by applying more strin- 
gent particle identification (PID) requirements to the bachelor kaon. To suppress back- 
grounds from B° — > D*~Df decays where — > K + 7i^7i + in the triple-bachelor decay 
B° D*~ K + 7t~tt + , it is required that m(K + ir~TT + ) is more than 15 MeV/c 2 away from 
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the -D+ mass. Reconstructed B° candidates are required to be well separated from the 
nearest PV, with decay time larger than 0.2 ps and good quality vertex fit. Candidates 
passing all selection requirements are refit with both D° mass and vertex constraints to 
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Selection efficiencies and trigger pass fractions defined below are evaluated using Monte 
Carlo simulation. In the simulation, pp collisions are generated using Pythia 6.4 [12] 
with a specific LHCb configuration 13 . Decays of hadronic particles are described 

The 



by EvtGen |14| in which final state radiation is generated using Photos 15 



interaction of the generated particles with the detector and its response are implemented 



using the Geant4 toolkit 16 as described in Ref. 17 



The simulated events are passed through an emulation of the hardware trigger and 
then through the full software trigger as run on data. The total kinematic efficiency, ekin, 
is determined from the simulation as the fraction of events that pass all reconstruction 
and selection requirements and the trigger. The fraction of selected events that pass the 
particular trigger selection relative to the total number of selected events is taken to be 
the trigger pass fraction, / trig . This fraction does not represent the true trigger efficiency, 
since it is evaluated with respect to a sample of events that have all passed the trigger. 



3 Fits to data 

The reconstructed invariant mass distributions for B° — > D*~h + 7r~7r + and B° — > D*~h + 
decays are shown in Fig. [TJ Simultaneous binned maximum likelihood fits to the CF and 
CS decays are performed, where the probability density functions (PDFs) are composed 
of a signal component and several background components. The CF signal shapes are 
required to share parameters with their CS counterparts. The total CF (CS) signal yield 
Acf (Acs) is given by the sum of the yield contained in the signal shape and the yield 
contained in the particle misidentification background from the CS (CF) sample. The 
values of Acf and the ratio Acs /Acf var y freely in the fit. 



The PDF of the signal decays is described by the sum of a Crystal Ball |18 and a 
Gaussian function, when all final state particles are assigned the correct mass hypothesis. 
This shape was chosen as it describes radiative loss and the non- Gaussian mass resolution. 
In the simultaneous fit to the CF and CS decays, the Crystal Ball widths vary freely and 
independently for both modes. A single freely varying peak position parameter is shared 
by all signal components in the fit to both modes. The Gaussian width is required to be 
equal to or greater than the Crystal Ball width, and this width is shared by the CF and 
CS signal shapes. 

All decay modes have background contributions from partially-reconstructed decays of 
the type B — > D*X , where X represents the final state bachelor(s) for the given decay 
plus an additional photon or pion that is not reconstructed. These backgrounds are 
parameterised by the sum of two single-sided Gaussian functions with a common mean 
and independent widths, all of which vary freely. A combinatorial background is present 
in all cases and is fit by an exponential function. The yields of the partially-reconstructed 
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and combinatorial backgrounds vary freely in all parts of the fit to data. 
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Figure 1: Invariant mass distributions of selected (a) B° — > D*~.frr + 7r~7r + , (b) 5° ->■ D*~K + , 
(c) -B —7- D*~7r + 7r~7r + and (d) I? — >• D*~it + candidates. Green (red) solid lines represent the 
CF (CS) signal shapes and their respective particle misidentification backgrounds. Orange dashed 
lines at low invariant mass represent backgrounds from partially reconstructed decays. Cyan 
dotted lines represent combinatorial backgrounds. The pink dash-dotted line below the signal 
peak in the fit to B° -> D*~ K + it 7r + candidates represents the background from misidentified 
B° -»• D*-Df and B° -> D*-K+K*° decays. 



Misidentification of pions and kaons causes cross-feed between the CF and CS signal 
decays. The simulation indicates that backgrounds of this type can be described by 
a PDF comprising two Crystal Ball functions that share a common mean but have 
independent widths. The absolute rate of a background from particle misidentification 
is fixed to be (100 — epm)% of the corresponding signal yield in the simultaneous fit, 
where epro represents the efficiency for all final state hadrons in the signal decay to be 
correctly identified. For the single bachelor decays, epro is determined by re-weighting the 
ALL(K — it) distributions obtained from calibration events to match the properties of 
the signal bachelor. The re-weighting is performed in bins of momentum, pseudorapidity 
and number of tracks in the event. Calibration tracks are taken from D* + — > D°7r+ decays 
with D° — > K~tt + [9]. For the triple-bachelor decays a similar approach is used, but the 
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Table 1: Selected candidate yields from fits to data that are used in the branching fraction 
calculations. The yield for a decay is given by the sum of the signal shape yield in the CF (CS) 
fit and the corresponding misidentification background yield in the CS (CF) fit. Uncertainties 
quoted are statistical only. 



Decay 


Yield 


B° D*- 


~7T + 7T 7T + 


7228 ± 93 


B° D*- 


-7T+ 


15,693 ± 136 


B° -> D*- 




519 ± 30 


B° -> D*- 


K+ 


1241 ± 53 



kinematic correlations between each bachelor are also considered. A per-candidate particle 
identification efficiency is determined from the product of each bachelor PID efficiency, 
and epiD is given by the weighted average of the per-candidate efficiencies. 

Further backgrounds are present in the B° — > D*~ K + -k~'k + sample from the decays 
B° —7- D*~D+ where Df — > K + K~tt + and the K~ meson is misidentified as a 7r~, 
and B° — > D*~K + K*° where K*° — > K~n + and the K~ is misidentified as a n~ . The 
backgrounds are modelled together using a Crystal Ball shape that peaks at a lower mass 
than the signal, with a peak position and width that vary freely in the fit. The Crystal 
Ball tail parameters are fixed to the values found in simulation. This background could 
be reduced by applying tighter particle identification requirements to the n~ bachelor, 
but this has not been applied in order to maintain symmetry between the CF and CS 
particle identification requirements. To eliminate the background from B° — > D*~ 
decays where — >• K + ti~'k + , the veto previously described is applied. The background 
from B° — > D*~Df decays with — > 7i + n^n + does not contribute due to the tight 
particle identification requirement applied to the bachelor kaon in the B° — > D*~ K + -k~'k + 
decay. 

The fits are superimposed on the data in Fig. [TJ and the measured yields Af for each 
decay are listed in Table [TJ The CF and CS ratios of branching fractions are obtained 
using 

M(B° D*-7i + 7i-7i + )/e tot (B° D*-7r + 7r-7r+) 
r 3h = fshX - 



Af(B° D*-n+)/e tot {B° D*-ti 

. . !\ 7V7T f l\ K7T , . ^ ' 

'cs — Jcs x 



*nr = f K** x ^(B° ~> D*-K+n-n+)/e tot (B° -» D^K+iTif 
cs Jcs Af(B° ^ D*-n+iT-iT+)/e to t(B ^ D*-it+it-it+ 

k rK M(B° -> D*-K+)/e tot (B° -> D*~K+) 

n cs — Jcs 



J\f(B° ->■ D*-7i+)/e tot (B ->■ D*-n+) ' 
The values of e to t are listed in Table [2j and the / factors correct for systematic effects. 
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Table 2: Kinematic efficiencies, trigger pass fractions and their product, taken from simulation. 
Quoted uncertainties come from the use of finite size samples to determine efficiencies and are 
accounted for as a source of systematic uncertainty. 



Decay 


Ckin (%) 


/trig (%) 


etot (%) 


B° -»■ D*- 


~7T + 7T 7T + 


0.037 ± 0.001 


69.3 ± 0.5 


0.0259 ± 0.0005 


B° -»■ D*- 


-7T + 


0.197 ± 0.002 


75.4 ± 0.3 


0.148 ± 0.002 


B° -»■ D*~ 




0.044 ± 0.001 


67.4 ± 0.9 


0.0298 ± 0.0007 


B° -»■ D*- 




0.201 ± 0.003 


75.4 ± 0.5 


0.151 ± 0.002 



4 Systematic uncertainties 

By measuring ratios of branching fractions, many sources of systematic uncertainty cancel. 
The systematic uncertainties for the r 3 h, r cs^ an< ^ r cs measurements will be discussed in 
turn. The primary sources of systematic uncertainty that remain in the r 3 h measurement 
are due to the different topologies of the signal and normalisation decays. Compared to the 
B° — > D*~7i + normalisation mode, the triple-pion decay mode has two additional pions 
which must be reconstructed and selected. The tracking efficiency has been studied using 
a tag-and- probe method with J/ip — )■ /i 4 ji~ decays [19] , which leads to a correction in r 3 h 
of 1.017 ± 0.035. In the B° — > D*~7t + 7t~tt + decay, the three bachelor pions are required 
to have a common vertex. The IP resolution and vertex x 2 /ndf distributions are observed 
to be ~15% broader in data relative to the simulation [I], resulting in a correction on r 3 h 
of 0.982 ± 0.016. 

Possible background from decays of the type B° — > D*~Df, where — > 7r + 7r~7r + , has 
been considered, and a correction of 0.990 ± 0.005 is applied to r 3h . The use of simulated 
events to determine trigger pass fractions has a residual systematic uncertainty arising from 
differences between data and simulation with respect to the emulation of the hardware 
trigger and trigger software. A correction of 1.009 ± 0.012 is applied to r 3 h to account 
for this difference. The candidate selection is limited to the mass region m(/T + 7r~7r + ) < 
3 GeV/c 2 . The m(7r + 7r~7r + ) distributions in data and simulation are in good agreement, 
such that the selection efficiency properly accounts for this choice. The fraction of events 
falling beyond 3 GeV/c 2 in simulation is 3.7%. Assuming 50% uncertainty on this value, a 
relative systematic uncertainty of 1.9% is assigned. 

The methods used to determine epio have an uncertainty from which the systematic 
contribution is determined to be 0.8%. A systematic uncertainty of 0.6% arises from 
the specific choice of PDF shapes in the fit. Both of the CF simulated samples have a 
comparable number of events after selection requirements are imposed, from which a 2.1% 
systematic uncertainty due to finite simulated samples is incurred. 

The CF and CS B° — > D*~h + 7T~7T + modes have identical selection requirements, apart 
from the particle identification requirements placed on the h + and the Df veto applied 
in the CS case. A systematic uncertainty of 3.4% is incurred as a result of the particle 
identification requirement placed on the bachelor kaon in the B° — > D*~ K + -k~'k + mode. 
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This tight requirement is necessary in order to reduce the background from misidentified 
B° — > D*~7t + 7t~7t + decays. To evaluate the loss of signal events due to the veto in the 
CS selection, the fit to data is performed both with and without the veto applied. The 
measured CS signal yield decreases by 1% upon application of the veto, which is taken as 
an inefficiency with 50% uncertainty and a correction of 1.010 ± 0.005 is applied to r^ n . 

The ratio of trigger pass fractions taken from simulation is 
/trig(#° -> D*-n + iT-iT + )/f tlig (B -> D*-K + it-it + ) = 1.03 ± 0.02, where the quoted 
uncertainty is derived from the size of the simulated samples. For such similar decay 
modes, this ratio should be close to unity. The ratio itself is already applied as part of the 
branching fraction calculation, but half of the difference in the ratio from unity (1.5%) is 
taken as an additional source of systematic uncertainty. 

In a similar fashion to the CF measurement, the CS measurement has a systematic 
uncertainty of 1% from the specific choice of PDF shapes and 3.0% uncertainty from the 
use of finite simulated samples to determine efficiencies. The fraction of CS decays with 
m(K + 7r~7T + ) in the range 2.7-3 GeV/c 2 in data is 5.5%. The fraction of events falling 
beyond the analysis cut at 3 GeV/c 2 is estimated to be half of this value (2.8%) with 50% 
uncertainty. The value of is therefore corrected by a factor 1.028 ± 0.023, where 
the quoted uncertainty contains a 1.9% contribution from the corresponding systematic 
uncertainty in the CF decay. 

The -B D*~h + modes have identical selection requirements apart from the particle 
identification requirements placed on the h + bachelor. A systematic uncertainty of 2.0% is 
incurred as a result of the particle identification requirements applied to the bachelors. 
Further systematic uncertainties of 0.7% , 1.7% and 2.0% arise from the specific choice of 
PDF shapes, the use of finite simulated samples and the trigger emulation 20 , respectively. 

Each contribution to the systematic uncertainty is listed in Table |3j The total systematic 
uncertainty is given by the sum in quadrature of all contributions. The overall systematic 
uncertainty for the CF measurement is 5.0%, with a factor f 3h = 0.998 that is applied 
as part of the calculation for r 3 /j. The CS triple- and single-bachelor measurements have 
overall systematic uncertainties of 5.4% and 3.4% respectively. A factor f^ 7 * = 1.038 is 
applied to r^g r7r , whereas Jq S = 1. 



5 Results 

The results for the ratios of branching fractions are 

r 3h = 2.64 ± 0.04 (stat.) ± 0.13 (syst.) , 
rg™ = (6.47 ± 0.37 (stat.) ± 0.35 (syst.)) x 10~ 2 , 
r£ s = (7.76 ± 0.34 (stat.) ± 0.26 (syst.)) x 10~ 2 , 

where the first uncertainty is statistical and the second is absolute systematic. Using 
the world average value for B(B° -»■ D*-tt + ) = (2.76 ± 0.13) x 10" 3 g, the branching 
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Table 3: Contributions to the relative systematic uncertainty for all measurements. The 
total uncertainty is obtained by adding the contibutions from the individual sources in 
quadrature. 





Value (%) 


Track reconstruction 


3.4 






Selection requirements 


1.6 






B° -> D*~D+ background 


0.5 






Trigger 


1.2 


1.5 


2.0 


m(/i+7r-7r+) > 3 GeV/c 2 


1.9 


2.2 




Particle identification 


0.8 


3.4 


2.0 


Choice of PDFs 


0.6 


1.0 


0.7 


Simulated sample size 


2.1 


3.0 


1.7 


Df -»■ i^ + 7T _ 7r + veto 




0.5 




Total 


5.0 


5.4 


3.4 



fractions are obtained 

B(B° -> L)*"7T + 7r"7r + ) = (7.27 ± 0.11 (stat.) ± 0.36 (syst.) ± 0.34 (norm.)) x 10~ 3 , 
B(B° -»■ D*-K + ) = (2.14 ± 0.09 (stat.) ± 0.08 (syst.) ± 0.10 (norm.)) x 10~ 4 , 

where the final uncertainty is due to the normalisation mode. Both re- 
sults are consistent with and improve upon the precision of the cur- 
rent world average values B(B° ->■ D*~n + n~n + ) = (7.0 ± 0.8) x 10~ 3 and 
B(B° -»■ D*-K + ) = (2.14 ± 0.16) xlO" 4 0. Combining the CF result and 
the current world average, where both values are weighted according to 
their total uncertainty, gives B{B° D*~tt + 'K~'k + ) = (7.19 ± 0.43) x 10~ 3 . 
The measurement of r^g 71 represents a first observation of the decay 
B° — > D*~ K + 7t~tt + . The value of this ratio is similar to the related measurement of 
B(B° ->■ D-K+it-it + )/B(B -> J D-7T+7r-7r + ) = (5.9 ± 1.1 (stat.) ± 0.5 (syst.)) x 10" 2 §. 
Using the updated world average value for B(B° — > D*~7t + 7t~tt + ), the branching fraction 
is obtained 

B(B° ->■ D*-K + 7T-7r + ) = (4.65 ± 0.26 (stat.) ± 0.25 (syst.) ± 0.28 (norm.)) x 10~ 4 . 



6 Search for excited charm resonances 

Using the same dataset, a search for orbital excitations of charm resonances (D**) con- 
tributing to the B° — > D*~7t + 7t~7t + final state is performed. The selection is identical 
to that presented above, except the lower px cut on bachelors is reduced to 300 MeV/c. 
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Figure 2: Corrected mass M(D*~ir + ) = m(D 
from background-subtracted B° — > D*~ir + ir~ 

contribution from 7Ji(2420)°, the red dashed lines represent contributions from the higher excited 
charm states D*(2460)°, £>(2550)°, 7J(2600) and 7J(2750), which are expected to be present [21 



but are not significant in this dataset. The blue dashed line represents B° decays that have not 
passed through an excited charm resonance. 



Events that have been triggered via tracks not associated with the B° — > D*~7t + 7t~tt + 
candidate are also included. The signal purity is only slightly reduced as a result of the 
looser selection. 

The corrected mass M(D*~n + ) = m(D*~n + ) — m(D*~) + m(7J*")p DG MeV/c 2 is com- 
puted for each D*~n + combination, in which the contribution to the mass resolution from 
the D*~ mass measurement is removed. To statistically subtract the background, each 



event is weighted using sWeights 22 obtained from the B° — > .D*~7r + 7r~7r + invariant mass 
fit. Fig. [2] shows the resulting distribution of M(D*~tt + ) for B° — > /J*~7r + 7r~7r + signal 
decays. 

A peaking structure associated to the 7J 1 (2420)° resonance is observed, where ~90% 
of the candidates in the peaking structure originate from the combination with the softer 
7T+ meson. Other resonances, consistent with 7J*(2460)°, 7J(2550)°, 7J>(2600) and L>(2750), 
are included in the fit but are not found to be significant. All resonances in the fit are 
described by Breit-Wigner functions. The 7Ji(2420)° Breit-Wigner function is convolved 
with a Gaussian resolution function of 3 MeV/c 2 width. The means and natural widths of 
all peaking structures vary around their established values [2l] with Gaussian constraints. 

The background from B° — > D*~7t + 7t~tt + decays that do not pass through an excited 
charm resonance is described by a function comprising the two-body phase space equation 
multiplied by an exponential acceptance function, e ~ aM ( D * 7r+ ). The shape parameter a 
and all yields vary freely. The branching fraction ratio is calculated by comparing the 
fitted 7J 1 (2420)° — > D*~n + yield with the total number of accepted D*~7t + 7t~7t + events 
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in the sample. A correction is taken from simulation to account for the acceptance, 
reconstruction and selection efficiency for events in the region close to m(.Di(2420) ) 
relative to the efficiency averaged across the full phase space. The ratio of efficiencies is 
/ = 0.91 ± 0.04. A systematic uncertainty of 10% is assigned to the choice of background 
PDF, which is determined by remeasuring the _Di(2420)° yield after shifting the fit range 
by ±50 MeV/c 2 . 

The measured yield is J\f(B° ->■ D 1 (2A20)°n + ir-) = 203 ± 42 and the total number 
of B° signal events after the looser selection is Af(B° ->■ D*-it + it-it + ) = 10,939 ± 105. 
Using these values, the ratio of branching fractions is obtained 

B(B° -)• (D 1 (2420)° ->• J D*-7r + )7r-7r+) . . . . 2 



B(B° ->• D*-n+n-n+) 



= (2.04 ± 0.42 (stat.) ± 0.22 (syst.)) x 10" 



where the numerator represents a product of the branching fractions 
B(B° ->■ J D 1 (2420)°7T-7r+) and B( J D 1 (2420)° ->■ D*~7r + ). The Wilk's theorem statis- 
tical significance of the _Di(2420)° peak is 5.9a, which becomes 5.3a when the systematic 
uncertainty is included. This constitutes the first observation of the colour-suppressed 
B° ->■ ^1(2420)°^+^- decay. 



7 Summary 

In conclusion, B° — y D*~h +r n~'n + decays have been studied using B° — y D*~h + decays 
for normalisation and verification. The branching fractions of B° — y D*~it + it~it + and 
B° — y D*~K + decays are measured, and the CS B° — y D*~ K + tt~it + and colour-suppressed 
B° — y Di(2420)°7r + 7r _ decays are observed. The final results are 

r 3h = 2.64 ± 0.04 (stat.) ± 0.13 (syst.) , 

= (6.47 ± 0.37 (stat.) ± 0.35 (syst.)) x 10" 2 , 

rg s = (7.76 ± 0.34 (stat.) ± 0.26 (syst.)) x 10" 2 , 

B(B° ->■ (/J 1 (2420)° -> D*-7t + )tt-7t + ) . . „ 

-i B(B° ^ D-n+n-n+) = ^ ± °^ ^ ± ^ X 10 " 

The results for and r^ s represent an improvement in precision, and the measure- 
ments of the decays B° ->■ D*~K+it-it + and B° ->■ (^1(2420)° ->■ L>*-7r+)7r-7r+ both 
constitute first observations. 
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